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We introduce a new class of models for interacting particles. Our construction is based on Jaco- 
bians for the radial coordinates on certain superspaces. The resulting models contain two parame- 
ters determining the strengths of the interactions. This extends and generalizes the models of the 
Calogero-Moser-Sutherland type for interacting particles in ordinary spaces. The latter ones are 
included in our models as special cases. Using results which we obtained previously for spherical 
functions in superspaces, we obtain various properties and some explicit forms for the solutions. 
We present physical interpretations. Our models involve two kinds of interacting particles. One of 
f--^ 1 the models can be viewed as describing interacting electrons in a lower and upper band of a one- 

dimensional semiconductor. Another model is quasi-two-dimensional. Two kinds of particles are 
^s^j , confined to two different spatial directions, the interaction contains dipole-dipole or tensor forces. 
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I. INTRODUCTION 



There is an intimate relation between group theory and certain one-dimensional exactly solvable systems Q, 0, 13, 
Q\ ' The radial part of the Laplace-Beltrami operator on symmetric spaces induces in a natural way an interacting 

one-dimensional many-body Hamiltonian with a characteristic gv~ 2 (x n — x m ) interaction between the particles at 
• positions x n and x m . Here, g is the coupling constant and the function v may be a sine, a hyperbolic sine or the 
identity, depending on the curvature of the symmetric space under consideration. These and similar systems have 
been studied first by Calogero and Sutherland HUH. They have much in common with the Brownian motion model 
studied by Dvson as early as in 1962 Other forms of the potential have been introduced, such as the Toda 

lattice or the Weierstrass function, which generalizes the original form of interaction. We refer to all models as 

Calogero Moser Sutherland (CMS) models irrespectively of the interaction potential and the underlying Lie algebra. 

The first proof of exact integrability of some CMS-Hamiltonians have been given in [l^. Later a more general 
proof has been given in 0,0] by very different arguments. In this context, we also refer to the work in Re f. flp . 

More recently, these models have been studied in the framework of supersymmetric quantum mechanics (171 Il8| . 
Although we work with supersymmetry as well, our approach is different from this. Generalizations to higher space 
dimensions 0,H3,0| have also been proposed. Extensive reviews are given in Refs. [22ll2^ |. 

Our supersymmetric construction extends and generalizes the group theoretical approach in ordinary spaces by 
exploring the relation between the radial part of Laplace operators on symmetric superspaces and certain Schrodinger 
operators: In some cases, i.e. for special values of the coupling constant g, the solution of the interacting particle 
■ Hamiltonian can be written as an integral over the classical matrix groups, the orthogonal, the unitary and the 
symplectic group. These groups are labeled by the Dyson index = 1, 2,4, respectively. The coupling constant q is 
a function of the Dyson index 0. Similar relations for Schrodinger operators exist also in superspace [24L |25l 1261 127| . 
A classification of matrix supergroups and more general of symmetric superspaces has been given in Ref. [28| . In this 
contribution, we introduce a labeling of symmetric superspaces in terms of a pair of numbers 02) akin to Dyson's 
index in ordinary space. This label may further be continued analytically in 0\ and 02 to arbitrary combinations 
[0i,02)- Our construction leads to a natural supersymmetric generalization of the CMS model for interacting particles. 
Hence, we arrive at a new class of many-body systems. They are likely to be exactly solvable in the allowed parameter 
region. 

Our construction goes considerably beyond the one by Sergeev and Veselov [2^, |3(], • These authors arrived 
at superanalogues of CMS models, starting from the underlying root spaces of the superalgebra. They also give a 
solution in terms of superanalogues of Jack polynomials. Their models however, depend only on one parameter and 
are therefore different from ours which crucially depend on two. Some of our models are related to the many species 
generalization of CMS models in Refs. |32t l33| . In contrast to our approach, the latter construction is ad hoc and it 
is not based on superspaces. 

The models we are investigating have been communicated in |34j . where emphasis was put on their interpretation 
and possible applications. Here we focus on mathematical aspects of the models. In particular the question of exact 
solvablity is discussed and exact solutions for certain parameters 0\, 02 are presented. 
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The paper is organized as follows: For the convenience of the reader we briefly compile some results for the models 
for interacting particles in ordinary space in Section \n\ Various supersymmetric generalizations of the models for 
interacting particles are presented in Section ITTll In Section llVl we find certain solutions by deriving a new recursion 
formula. In Section we give an extensive interpretation of the physical systems described by the supersymmetric 
models. A brief version of this section can be found in [34(. We summarize and conclude in Section IVT1 

II. MODELS FOR INTERACTING PARTICLES IN ORDINARY SPACE 

In Section Hi Al we sketch the connection between ordinary groups and the many-particle Hamiltonian. We discuss 
the connection to the recursion formula in Section III Bl 



A. Differential Equation and its Interpretation 

The connection between some models of the CMS type in ordinary space and some radial Laplaceans appearing in 
group theory is seen by considering the eigenvalue equation 

A^f( X ,k) = -^2f^N(^k). (1) 

The N variables x n , n = 1, . . . , N are interpreted as the positions of the particles later on. There is a further set of N 
variables k n , n — 1, . . . , N which will play the role of quantum numbers. The operator depends on a parameter 
(3 and is given by 

n— 1 1 

where 

A N (x) = JJ (x n - x m ) (3) 

is the Vandermonde determinant. If the symmetry condition <$>^\x,k) — (fc, x) and the initial condition 

(0, k) — 1 are required, the solution of the eigenvalue equation is for /3 = 1,2,4 equivalent to group inte- 
grals over 0(N), U(iV) and USp(27V), respectively. These integrals are referred to as spherical functions 35]. We 
notice that they are different from the group integral which Harish-Chandra investigated in Ref. El Q . This is re- 
flected in the operator which is the radial Laplacean on symmetric spaces with zero curvature [36], more precisely 
on the spaces of symmetric, Hermitean, and Hermitean selfdual matrices for [3 — 1,2,4. Only for /3 = 2, the Laplacean 
coincides with the Laplacean over the algebra of the group U(iV). This is the only case where the spherical function is 
identical to a Harish-Chandra group integral due to the vector space isomorphism of Hermitean and anti Hermitean 
matrices. For arbitrary f3 the eigenvalue equation Q is closely connected to models of one dimensional interacting 
particles. Using the ansatz 

Ni ' k) -^\x)^(k) (4) 
the eigenvalue equation Q is reduced to a Schrodinger equation 

^(f- 1 ) E fa -x )2 WW) ( 5 ) 

^ ^ n<m ^ n m ' J \n— 1 / 

which contains a kinetic part and a distance dependent interaction. Often, one adds N confining potentials to the 
interaction in Eq. 0. This is done to make the system a bound state problem. However, apart from this, the structure 
of the model is not significantly affected by this modification. Thus, we will not work with confining potentials in the 
sequel. The specific model Eq. (0 is also called rational CMS model or free CMS model. 




3 



The solution 'Jjy (x, k) is now interpreted as a wave function of the Schrodinger equation JSJ) with energy ^ k„. 
Thus, no symmetry condition such as ty^p (x, k) = ^^(fc, x) is imposed. In the following we always refer to functions 
such as ^^\x,k) as wave function. On the other hand, functions such as &^\k,x) and more general solutions of 
eigenvalue equations of type are referred to as matrix Bessel functions. 

The parameter (3 > measures the strength of the inverse quadratic interaction. The interaction can be attractive 
(3 < 2 or repulsive f3 > 2. For (3 = 2, the model is interaction free. This is group theoretically the unitary case and 
equivalent to the Itzykson-Zuber derivation [37j of the U(iV) Harish-Chandra integral. 

The symmetric spaces mentioned above stem from a common larger group, namely the special linear group. In 
Cartans classification they are referred to as A, AI and All 5]. There are other symmetric spaces derived from the 
orthogonal and the symplectic groups as larger groups, designated B, C and D, respectively. These symmetric spaces 
are also related to Schrodinger equations, but with a different interaction |22^ . 



B. Connection to the Recursion Formula for Radial Functions 

For arbitrary positive (3 the solutions of the eigenvalue equation <£Q) $^' ) (x, k) can be expressed in terms of a 
recursion formula j^S Ei3 

^\x,k) = J d^(x',x) exp (j>„ - ®N-i( x ',k) , (6) 

where k) is the solution of the Laplace equation ||TJ for N — 1. Here, k denotes the set of quantum numbers 

k n , ri = 1, ... , (N — 1) and x' the set of integration variables x' n , n = 1, . . . , (N — 1). The integration measure is 

dn(x',x) = of An -_\ { ^ -n^-o d ^ ■ w 

A% (x) \ 11 ' 



a .rn 



Here, d[x'] is the product of all differentials dx' n , n = l,...,(N — 1). The constant G^ guarantees a proper 
normalization. The inequalities 

x n < x' n < x n+ i , n= 1, ...,(N- 1) (8) 

define the domain of integration. An equivalent recursion formula exists also for the eigenfunctions (x, k) of 
the Hamiltonian in Eq. (J5J. For f3 = 1,2,4 the above recursion formula is equivalent to group integrals over O(N), 
V(N) and USp(2iV), respectively. The case of arbitrary (3 has not found a clear group theoretical or geometrical 
interpretation yet. However, many properties which are obvious for the group integral carry over to arbitrary f3. 
We just mention the following. $^ (x, k) is a symmetric function in both sets of arguments. This has as a direct 
consequence that the behavior under particle exchange of the wave function ^E'{x, k) is only governed by the 

/2 3 1 2 

Vandermonde determinant A^; (x)A^y (fc). The wave function obtains under particle exchange a complex phase 

Pum^Hx, k) = exp(-inl3/2)^{x, k) , (9) 

with 

P-nm^ w (^1) ■ ■ • i ^n? ■ • • j . . . , fc) — ^ N (^1' ■ • * ? ^mi • • • : ^ni . . . , fc) . (10) 



For this reason the model of Eq. JSJl is frequently used as paradigm for systems with anionic statistics [4(J, |4l| . A 
recursion formula akin to Eq. (jSJ has also been derived for Jack polynomials |42|. 



III. MODELS FOR INTERACTING PARTICLES IN SUPERSPACE 



A classification of supergroups and superalgebras similar to Cartan's classification in ordinary space can be found in 
Refs. |43l l44j . Apart from some exotic groups, there are essentially only two families of supergroups. The general linear 
supergroup GL(fci/fc2) respectively its compact version the unitary supergroup U(fci /ki) and the orthosymplectic 
group OSp(fci/2fc2). A classification of the symmetric superspaces has been given in Ref. 28]. 
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In Sections IIII Al and IIII Bl we present supersymmetric generalizations of models for interacting particles based on 
the supergroups GL(fci/fc 2 ) and on the symmetric superspaces GL(/ci/2fc 2 )/OSp(fci/2fc 2 ). In Section UlI PI we give 
the supersymmetric generalization based on the supergroup OSp(fci/2fe2)- In Sections IIII CI and IIII El we introduce 
two more general models which comprise the other models derived before as special cases. These models can be 
considered as supersymmetric generalization of the Schrodinger equation J5| for the CMS models in ordinary space. 



A. Models Derived from the Superspace GL(fci/fa) 

To extend the models in ordinary space to superspace, we begin with models derived from the superunitary case. 
The underlying symmetric superspace is called A|A in Ref. |2jj|. We construct the eigenvalue equation 



A^Ag 2 ( S) r) = ^ + £ r 2 2 ) xflM . (11) 



for the operator 

Z,. 1~ _ 

d 



U B U (-) klk2 B * VP U < fc2 (-) 9s p2 ^ 



"* w *£ • (12) 



where the function I24LI4 



, , _ IIp< 9 (^l - Sql) U p<q (Sp2 - S q2 ) 

B kl kA s ) - ff — / : \ 



is the square root of the Berezinian for the superalgebra u(fci/fc 2 ). Using the ansatz 
leads to the Schrodinger equation 



A *ifal«>»V - p/2 0/2 (14) 



\p=i P =i / 

which includes the eigenvalue equation as special case for k\ = or fc 2 = 0. Again, the case (3 — 2 gives, for 
all k\ and fc 2 , an interaction free model, connecting to the supersymmetric Harish-Chandra integral for the unitary 
supergroup U(fci/fc 2 ). 



B. Models Derived from the Symmetric Superspaces GL(fci/2fc2)/OSp(fci/2A;2) 

Also the two forms of the symmetric superspace GL(fc!/2fc 2 )/OSp(fci/2fc 2 ) yield new supersymmetric models as 

well. These spaces are denoted AI|AII and AII|AI in Ref. [2^|. They involve the Berezinians B^ 2k2 (s), see Ref. [25$ . 
Apart from some absolute value signs which are not important here, one has c = +i for the symmetric superspace 
AI|AII and 



^ rip< 9 (^i-^i)n P < 9 (^2-s g 2) 4 
Bk ^ {s) = rw^-M 2 (16) 

while one has c = — i for the symmetric superspace AII|AI and 

S(-l) ( , _ U p<q (s P l-S ql )U P<q ^ P 2-S q2 ) 4 
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Thus, we obtain the radial part of the Laplace-Beltrami operator 

k 2 



A (c) _ \ - 1 9 s(°) ^ d l r 1 9 5 (c) , , d 

s ~ h w 5 ^ fei2fc2 ( } 9 v 5 p ^ ( S) ^ 2 Bfci2fe (s) a* (18) 

and the eigenvalue equation corresponding to Eq. 



Employing the ansatz 



2 

Vp=l p— 1 



<u^)= ~ (e) . Ki ;~ c ' , (20) 



we find the Schrodinger equation 



k\ ~n . /c 



^72~ + oL TfTa - + o 7 \2 ~~ 2 7 72 ~ 7 72 I "fcifc2 ( s ' r ) 



E^+Efl 41 m- (2D 

\p=l p =i / 



The choices ki = and fci = in Eq. I|21ll yield Eq. (JSJ with /? = 4 and j3 = 1, respectively. For arbitrary ki 
and ^2 the function (s, r) is the supersymmetric generalization of spherical functions which we treated in a 



previous work |26l l27fl. For fci/2 = k2 — k these models are of prominent interest in random matrix theory. The 
fc-point eigenvalue correlation functions for a random matrix ensemble can be expressed as derivatives of a generating 
functional. This generating functional obeys a diffusion equation in supermatrix space |25j similar to Dyson's Brownian 
motion in ordinary matrix space 0, ITo| . The kernel of this diffusion equation is given by the solution of Eq. (|21|) . 



C. Embedding of the GLQti/kz) Based Models into a Larger Class of Operators 

We now embed the functions Bk 1 k 2 { s ) an d B^k ( s ) °f Eqs. <|13[) . (|16l) and (|17() into a larger class of functions 
defined by 



n Pig (%i-c Sg2 )^ 



4X 2J ( fi ) = llp<qK _ — P<q 7^ — ■ (22) 



Here, we introduce two parameters fa and fa. This is of crucial importance for the resulting models. They become 
very rich due to this twofold dependence. We assume that these parameters are positive, fa, fa > 0. The parameter 
c can take the values c = ±i. The functions B^'^ 1 '^ 2 ' (s) induce a differential operator 

A (cA,/32) _ 1 V 1 9 B {c ' PuP2) (c) 9 I 1 V 1 8 B (c,fa,fa)/ v g f23 s 

V/^ a^A*- (s W + v/^ s£f >(«) ^ felfe2 ( W ■ ( j 

In the first quadrant of the (fa, fa) plane B^ > '^'^ 2 \s) and therefore Ai c ' /3l ^ 2 ' ) is analytic in fa and fa, respectively. 
The eigenvalue equation corresponding to Eq. (fTT)l reads 

A^«it A) (.,r) = ~ (E ^ + E t|=) "fcfW) • ( 24 ) 



With the ansatz 
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we obtain the Schrodinger equation 




In the sequel, we refer to the model (|26|) as superunitary model. 

The superunitary model includes the models derived from the unitary supergroup, discussed in Section IlII Al for 
01 = /?2 = 0- The models derived from the symmetric spaces AI|AII and AII|AI discussed in Section IIII Dl arc 
included. They result for 0i = 1, 02 = 4 in the case c = i and for 0i = 4, 02 = 1 m the case c = — i. The solutions 

Pk^' M and ^i C 1 fc2 1 ' /32) are real anal y tic functions in ft and 2 . Since Ai - ^ 1 '* = A^ 1 '^ the solutions also have 
the symmetry 

*££ A) *(-i.-*.') =Cf lA) (»i,-a,r) = <£'*>(.i,-*,r) . (27) 

We observe that only in the case 0i = 02 = the interaction between the two sets of variables vanishes. If we choose 
ft = and /?2 7^ 0, we recover the noninteracting model, i.e. the Harish-Chandra integral, for the variables r pl , s pl , 
p = 1 . . . kx. Analogously, the choice ft = and ft 7^ yields the noninteracting model, i.e. the Harish-Chandra 
integral, for the variables r p2 ,s P 2, p = 1 . . . k 2 . In Eq. (|2^|l the points {fa, fa) = (0,0) and (01,02) = (2,2) are 
indistinguishable. They both yield a completely noninteracting model in either set of variables. As mentioned before, 
the point (2, 2) has the group theoretical interpretation as supersymmetric Harish-Chandra integral. 

The CMS models in ordinary space Eq. (JSJ are recovered by setting either k\ = or k 2 — 0. For the models Eq. (JjjJ 

the points of even — 2,4, 6, . . . are special )l6j . The wavefunction can always be written in an asymptotic 
expansion akin to the Hankel expansion of Bessel functions |46|. In a previous publication |38ll39| we showed that only 
for even this asymptotic expansion terminates after a finite number of terms. In the present context, this property 
carries over to the points (0i,0 2 ) = (2n, 2n), n £ N+, since there the Schrodinger equation Eq. (|26|l decouples into the 
sum of two independent CMS models Eq. J5| . It is an intriguing and unsolved question if there are other points in the 
(0i, 02) plane with this property. We conjecture that this property holds for an arbitrary point (2n, 2m), n, m € N + . 

Due to the non-Hermitecity of the left hand side, the interpretation of Eq. i|26|) as a Schrodinger equation has to 
be done with some care, see Section Ivl 



D. Models Derived from the Superspace OSp(fci/2fc2) 

Furthermore, we derive another class of models by considering the group OSp(fci/2fc2) instead of GL(fci/fc2). The 
role of the Berezinian Bf, 1 2k 2 ( s ) is taken over by the functions |47J 

CkMs) - n„>2i+^) (8) 



for even ki and by 



n i \ n p<g ( S pl IIp<g( S p2 ^2) Up=l 1 S pl , , 

CkMs) = n M (-2i+^) (29) 

for odd ki. Here, we employ the notation [fci/2] for the integer part of ki/2. The two formulae differ only in the last 
terms of the numerators. We define the operator 



A(uosp,/3) _ 1 1 d p f ] 9 , ] f 1 d C P l s \JL_ 

2 t^l C k,2k 2 ( S ) dSpl dSpl P=l C kl 2k 2 ( S ) dS P 2 dS P 2 



(30) 
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such that we recover the supersymmetric Harish-Chandra case for (3 — 2, see Ref. [42|. We seek the eigenfunctions 
X*fL 2 ( s ' r ) of this operator, 



%/2] k 2 

p— i p— i 



To arrive at a Schrodinger equation, we make the ansatz 



which yields 



,09) 



Cfi / 2fc 2 ( S ) ( ^fei / 2fe2 ( r ) 



,/3/2 



(31) 



(32) 



2 £s »« 



-. fc 2 



2 



fci2/c 2 1 



[fei/2],fe 2 




a pl,2 



(33) 



The last sum on the left hand side of Eq. Q33|l extends over the variables s P 2,P = 1 . . . Afc in case of the Berezinian <|28l) 
and over s p i,p — 1 . . . [fei/2] in case of the Berezinian l|29[l . Once more, we arrive at an interaction free model for 
(3 = 2, corresponding to the supersymmetric Harish-Chandra integral over the supermanifold UOSp(fci/2£;2), see 
Ref. Ui. A gain, as before in case of the unitary supergroup, there is no interaction between the two sets of variables 
s p i and s P 2- This is so for all values of (3. We notice that for arbitrary (3, the model introduced here contains two 
models in ordinary space which are not included in the models of Section [H] For ki = 0, we obtain models based on 
O(ki) and for k\ = 0, we obtain models based on USp(2/c2)- Both were discussed in detail in Ref. |2^ . 



E. Embedding of the OSp(fci/2fc2) Based Models into a Larger Class of Operators 



In Section llll CI we embedded the models of Sections llII Al and lHI Bl into a much richer structure with two parameters 
Pi and /?2- We now perform the analogous embedding for the OSp(fci/2fc 2 ) based models l|3^|) . Here, the result is 



[fcx/2] 



-y 



d 2 



pi 



L v — 



(h 

2 



E 



24. 



p2 



2^2 



2i 

2 



E 



Hi + Hi 



[ki/2 



p<q ( S pl S ql) 



(1-0E 
1 



2 s 2 
n=l /S « 2 



>fil 



{ 2<r 



;V / A& + l)E 



'pi 



P,9 ( S pl + S 92) 



p,q S Pl 



S «2 / 



/[fcx/2] 1 te 1 \ 



(s,r) 



(34) 



We introduced the quantity Z with / = for k\ even and I = 1 for k\ odd. In the sequel, we refer to the model 1341) 
as orthosymplectic model. 



For (3x = (h = (3, Eq. I|33J) is recovered from the orthosymplectic model with 



^l2fe 2 



(s,r) 



(s,2r). The 



discussion of Eq. (|34|l is along the same lines as the one at the end of Section IIII CI For k\ even — in analogy to 
the model based on the unitary supergroup — the points (fix, (h) = (1,4) and (/3i , /?2 ) = (4, 1) correspond to certain 
symmetric superspaces, namely to the two different forms of the symmetric superspace OSp(fci /2fc2)/GL((fci /2) /k^). 
They contain the symmetric spaces SO(fei)/SL(fci/2) and Sp(2fc2)/SL(fc2) as submanifolds. In Ref. [23 they are 
denoted CI | Dili and DIII|CI, respectively. 
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IV. SOME SPECIFIC SOLUTIONS 



The superunitary model (|26|l and the orthosymplectic model i|34|) . comprising the GL(/ci//c2) and the OSp(fci/2fc2) 
based models, respectively, have a very rich structure due to the dependence on the two parameters fa and fa. Thus, 
the general solutions are highly non-trivial and not known to us at present. Nevertheless, we are able to construct 
exact solutions of the superunitary model given in Eqs. l(2"4"|) and (|2"5|l for special values of the two parameters {fa, fa). 
More precisely, we derive solutions on certain one-parameter subspaces of the (fa, fa) plane. We distinguish two such 
one-parameter subspaces: first, the diagonal fa = fa and, second, the hyperbola fa = 4//?i, see Fig.^ The solutions 
in these subspaces contain the solutions of the models introduced in Sections IIII Al and IIII Bl 




FIG. 1: Curves in the (fa, fa) plane for which we construct solutions. The group theoretical points are also indicated. 

In Sections ITVAl and ITVBl we state and discuss the solutions on the diagonal and on the hyperbola, respectively. 
The solutions we derive on the hyperbola are generalizations of the recursion formula stated in Section III Bl In 
Section TlV CI we give the derivation. We also present some few-particle solutions in Section TlV Dl 



A. Solutions on the Diagonal /3 = fa = fa 

In this case the Schrodinger equation is Eq. (|15|l . It decouples into equations for two independent sets of particles. 
Hence we can write the solution as the product 

vi% 2 (*,r) =*g ) ( S i,n) ■ (35) 

An explicit expression of (sj, r i ) 1 i = 1,2 can be obtain in terms of the recursion formula Eq. © and in 
combination with Eq. We point out that the calculation of each ^^\si, separately is in itself very difficult. 



B. Solutions on the Hyperbola /3 = fa = 4/ 'fa 

In general, the interaction term between the two different sets of one-dimensional particles at positions s pl and s q 2 in 
the superunitary model Ij26(l does not vanish. Thus, as the product form of Eq. 1|35JI is destroyed, it is highly non-trivial 
to obtain solutions for positive parameters Pi and fa and for arbitrary dimensions k\ and &2- Nevertheless, we derive 
solutions on the hyperbola fa = 4//3i shown in Fig. ^ On this hyperbola, the exponent in the denominator of the 
function (1221) has the constant value two. Exactly on the same hyperbola fa = 4/ fa, Sergeev and Veselov constructed 
supersymmetric extensions of CMS models and their solutions in terms of deformed Jack polynomials |3(l l31| . The 
existence of a recursion formula suggests that there should be a recursion formula akin to the formula derived by 
Okounkov and Olshanski f° r th e deformed Jack polynomials as well. 

We emphasize again that we expect the superunitary model l|26|) to be exactly solvable for any positive fa,fa- 
As argued in Refs. |38ll39||. the recursion formulae can be viewed as generating functions for Jack polynomials, 
or equivalently, as a proper resummation. This carries over to the present case. We generalize the supersymmetric 
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recursion formula of Ref. for the symmetric superspace AI| All discussed in Section llll Dl We analytically continue 
the solution at the point (f3i,fl 2 ) = (1)4), (4, 1) and (2,2). Thereby we construct the solution on the hyperbola. 

We write p^'jf = Pkits ^ or ^ ne s °hition of the superunitary model (|24|l on the hyperbola. In Section TlV CI it 
will be proved that p*£'fP can be expressed through the recursion formula 



fei fel — 1 

\p— l p— l p— l 

Here, pfe^i)fc 2 (s',?*) is the solution of the superunitary model (|24|l for the fei + k 2 — 1 variables s' = 
[s'i\, ■ ■ ■ , s '(fc 1 _i)i, s 'i2) ■ • • ) s fe,2)- The solution is labeled by the quantum numbers r — (rn, . . . ,r(k 1 -i)i}i'i2, ■ ■ ■ , r fc 2 2)- 
The primed variables are the integration variables. The integration variables s p i are commuting. Their domain of 
integration is compact and given by 

s P i < s' pl < S(p+i)i , p = 1, . . . , (fei - 1) . (37) 

The integration variables s' p2 are related to Grassmann variables £ p and £* by 

ICpI 2 = csp 2 - c Sp2 ■ (38) 
The modulus squared of a Grassmann variable is defined by 

I£pI 2 = C P tv = - ZpCp , (39) 

which is the formal analogue for the length squared of a commuting variable. The integration over Grassmann variables 
is defined by 

[d£pdt; = and J \i P \ 2 di p dC p = 1 • (40) 

The normalization to one differs from the convention we used in Ref. |2?|. where the integral was normalized to 1/2tt. 
The integration measure d^ c '^(s', s) reads 

^(s',s) = ^ ) (s' 1 ,s 1 )^' f3 \s' 2 ,s 2 )^f(s',s) , (41) 
with the products of the differentials 

k 2 fci-1 

= n <%p<%* p and d i s i] = n Ki » ( 42 ) 

p— i p— i 

and the measure functions 

/ \ 2/0-1 

,WfJ „ \ _ A /„'U 1 — 4/ /3 / 



<Vi,*i) = A fcl ( s ;)A^( Sl ) -II 



p.'-; 



^(4*2) = AC /4 (- 2 )AC /4 ^(- 2 )fi(^- c 4 2 )" /2 ^ 2/4 

v>bf\ s> > s ) = nnn^-^v-^^-v^-^ ■ (43) 

p=i ;=i g=i 

We split the measure function into three parts /l«b, ^_f, ^bf as in Ref. [27| • We do so, because the coordinates are 
originally for certain values of 0i and [3 2 , Bosonic and Fermionic eigenvalues of some supermatrices. The recursion 
formula (|36(l reproduces the recursion formula derived in Refs. [2"fiLl27j | for /3 = 4. It also reproduces the supersymmetric 
Harish-Chandra integral discussed in Ref. |4j| for (3 — 2. Moreover, for k 2 = Q the recursion formula in ordinary 
space found in Refs. |38l l39j| and briefly discussed in Section III Bl is naturally recovered. 
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The case ki = deserves some special attention, because [Ib and \ibf vanish and so does the exponential in Eq. I|36|l . 
Importantly, the function /ip does not. The corresponding Schrodinger equation is just that of the CMS-Hamiltonian 
for fc 2 particles as defined in Eq. ®. Its solution, or more precisely the solution of its associated Laplace equation 
is by definition given by Pgfca = • However, the recursion formula yields another solution 

Sg> (s 2 ,r 2 ) = J d[^ F ( S2 ,s' 2 )^( S ' 2 ,r 2 ) . (44) 

For this to hold the Laplacean A if ' defined in Eq. @ has to commute with the Grassmann integration of Eq. 14411 . 
This implies that the eigenvalues of the operator defined through the Grassmann integration are conserved quantities. 
Indeed, the operator Aif 1 commutes with the Grassmann integral Eq. H44JI . 

Aif J d[^ F (s2,s' 2 )f(s' 2 ) = I d[^ F {s2,s' 2 )^f{s' 2 ) , (45) 

where f(s' 2 ) is analytic and symmetric in its arguments, but otherwise an arbitrary test function. We sketch the 
derivation of Eq. lf4l)|l in E] 



C. Proof of the Recursion Formula 



We now prove that the functions p^kl gi ven by Eq. indeed solve the differential equation Eq. (|23|l on the 
hyperbola. The proof relies on the invariance properties of the measure function fj,( c '^(s' , s). We define the Laplace 
operator Eq. (|23[) on the hyperbola Ai c '^ = /sjf' 4 ^ 13 ' 13 ^ anc j the center of mass momentum operator 



fci 



pi e) - E 



, dspi 



&2 n 



(46) 



We then have the two identities 



>(c) 



d/J,(s,s')f(s[,s' 2 ) = I d/i(s,s')P s ( , C) /(si ; 4) 



A^J M*,*')M,s'*) = J d^s lS ')A^f( s [,s' 2 ) , (47) 

which hold for an arbitrary function f(s\, s 2 ) symmetric in both sets of arguments s p \, p — 1 . . . ki and s P 2, p = 
1 . . .k 2 . We derive Eqs. (|47|l by direct calculation, using repeated integration by part. This procedure is relatively 
simple for the first equation of l|47|l . However, for the second one it becomes rather tedious due to the complexity of 
the measure function. Some of the steps are sketched in|5| A more elegant proof is likely to exist. 

Employing the properties l|47(l , we can now prove the recursion formula by acting from the left with Af' on both 
sides of Eq. We set 



/(si, 4) = exp (-i s ' P i ~ f E* s ^ rfe i^ ^(fcf-i)fe 2 ( s '' ? ) 



(48) 



and obtain straightforwardly from 147|l 



fcl-l 



E S P! ~ E S Pl + f E l^l 2 ) rfe l! 



0=1 



p=l 



p=l 



-^-rkil+^s' I P(k 1 -l)k 2 ^ ' ' 



Since by definition we have 



A (c,/3) (c,/3) , , ~ 



E 



fc 2 „2 



Y pl 

=i p=i 



E^)4V(^ 



(49) 



(50) 



11 



we arrive at 



^& 2 (s, r) = - ( £ ^ + e % ) « (s , r j , ( 5D 
\ p =i 



2 ' ^ ^ 

p—i v ' / 



which is our assertion. 



D. Few Particle Solutions 



Once the eigenfunction ^(s2, r 2 ) in ordinary space is known, we can recursively construct the eigenfunctions 



?( s ' r ) lrom formula l|3l)|l by starting with p ( Q^'{s, r) = <& ( ^ p> {s 2 , r 2 ). The eigenfunctions <J>j£ p,l (s2, r 2 ) are given by 



(c,/3), 



the recursion formula in ordinary space, see Section [II Bl We illustrate the procedure for two examples in superspace. 
For the sake of simplicity, we consider only c = +i and suppress the upper index (c) in the sequel. 
To begin with, we study the case k\ = k 2 = 1- The eigenvalue equation is 



1 £P_ 



yielding the closed solution 



p^' M (s,r) = exp 



1 (P_ 



i 



Sll - «Sl2 







ds n 



- d 

''ds 2 i 



' 11 
/3I 



' 12 



(52) 



±- 



(\/A" s ii - i\/Jhsi2^ (r n - ir 12 ) 



Iv^-v^I^'VhJO*). (53) 



fpx-y/P2 

Here, H J/ (^) is the Hankel function of order = •^//3 1 /3 2 /4 + 1/2. Its argument is the dimensionless complex variable 



Oh 



(si 



IS12) 



(54) 



The result (|53[1 holds for all arbitrary positive parameters (3\ and (3 2 . 

From Eq. (|53|l we can gain deeper insight into the structure of the solutions on the hyperbola j3ij3 2 = 4. The order 
f of the Hankel function becomes 3/2 on the hyperbola. The asymptotic Hankel series of the half integer Hankel 
function of order n + 1/2 terminates after the n-th step ^(|. On the other hand, the asymptotic series of a Hankel 
function whose order is not half-integer is infinite. Thus, only the Hankel functions of half-integer order can be 
expressed as a product of a finite polynomial and an exponential. The value v = 1/2 corresponds to either f3\ = or 
fh = and hence to a one-type-of-particle model, see Eq. Q and Eq. JSJ. Consequently, the order v = 3/2 is the 
lowest half integer order describing a two-type-particle model that has a non-trivial solution which can be written as 
product of a polynomial and an exponential. Furthermore, we notice that it is exactly this extra term in the Hankel 
expansion of Hj, 2 (z) which can be expressed by an integration over properly chosen Grassmann variables. Indeed 
the recursion formula yields directly 



Pil (sn,si2,r 11 ,(3r 12 /2) = exp(ir n su + i/3r 12 si 2 /2) 



i-1 
2 



i(3 r 



sn) (irit - rn) 



(55) 



which is identical to Eq. I|53|l on the hyperbola. We expect recursive solutions of the k\ + k 2 particle Hamiltonian 
Eq. 1|26|) akin to the recursion formula Eq. I|36|) to exist for other half-integer v as well. 

The next simplest case is k\ = 1 and k 2 = 2 and vice versa. It is still possible although cumbersome to find an exact 
solution for arbitrary f3\ and [3 2 . As we only wish to illustrate how the recursion works, we do not derive this exact 
solution here. Rather we use formula (|36(l to find a solution on the hyperbola. Without loss of generality we choose 
k 2 = 2 and k\ = 1. The bosonic measure pb(si, s[) vanishes. We have only to perform four Grassmann integrations. 
This implies that the solution can be written as a differential operator acting on p£ 2 ^ {s 2 , (3r 2 /2) = (s 2l (3r 2 /2) 



, (s 1 ,s 2 ,n,l3r 2 /2) = LW(a, r)$f (a 3> /3r a /2) 



(56) 
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Using the definitions of the measure Eq. (|41() and Eq. (|43|l and doing the Grassmann integrations we find 



2 r 



n (is P 2 - s u ) < yi 



.9=1 



- 1 



1 (3 d 

^«77 r ll 

lSq2 ~ Sn 2 OS q2 



(3 f(3 







1 



2 V 2 / ]lp=l (« S P2 - Su) 2 (Si2 - S 22 ) V 5s 12 9s n 



(57) 



For the eigenfunction $2 ( s 2i l3r 2 /2), we employ the explicit form 



*(/3)/ a /0 \ / - fl ( s 12 +S22)(?-12 +r 22 )\ (B+i) { @z\ 

Q (s2,Pr 2 /2) = exp I -tj3 - I X u ' I -j- I 



which involves the spherical functions 



x (^l) H=2 (^D/2 r((/3+1)/2) ^g 



(/3-l)/2 



(58) 



(59) 



where J„ is the Bessel function of order v. The variable z — (si2 — S22)(?"12 — r 22) in Eq. (|58() is dimcnsionless. 
Plugging this expression into Eq. I|56(l and using Eq. 157JI we can cast pf^ into the form 

Pi2 (su, s 2 ,rn, (3r 2 /2) = exp (irusu + —i (ria + r22) (su + S22) 



^-1 
2 

/? 2 



^ 12 ^22 



/3-l + 2^+i^(ru-^-^] ( Sll -^-^ 



«Si2 «S 22 



■ n( r " - irjaHsu - is p2 ) X (f3+1) 
p=i J 



(60) 



which explicitly shows the symmetry between the two sets of arguments s and r. 



V. PHYSICAL INTERPRETATION 



To develop an intuition for the physics of the differential operators (|26|l and 13411 in superspace, we recall the 
physical interpretation of CMS models in ordinary space. The Schrodinger equation models a system of N 
interacting particles in one dimension, moving on the x-axis, say. The eigenfunctions are labeled by a set of conserved 
quantities or, equivalently, quantum numbers k n , n — 1 . . . N . This is tantamount to saying that the system is exactly 
solvable. In the limit of vanishing coupling, i.e. for (3 = 2, the quantum numbers are the momenta of each particle. 
The characteristic feature of this model is the (x n — x m )~ 2 interaction potential. The models based on the ordinary 
groups O(N) and Sp(27V) fit into the same picture. However, the models have in this case a symmetry under point 
reflections about x — - 0. Moreover, for the symplectic group and the orthogonal group with N odd, there is an 
additional inverse quadratic confining or deconfining central potential ■ 

We now show that the physical interpretation along those lines carries over to our superspace models in a most 
natural way. We discuss the superunitary model in Section IV Al and the orthosymplectic model in Section IV Bl 



A. Superunitary Model 



The superunitary model is given by Eq. I|26|l . We notice that its differential operator is not Hermitean. This leads 
to some ambiguity in the interpretation of the model. The imaginary unit in the parameter c is due to a Wick-type-of 
rotation of the variables s p2 . This was needed in Ref. |48j to ensure convergence of integrals over certain supermatrices. 
However, in our application, there is no such convergence problem, as long as we do not go into a thermodynamical 
discussion of the model. Thus, we undo the Wick rotation by the substitution is p2 — > s p2l p = 1 . . . k 2 . We introduce 
the coupling constants 
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.922 = V 02 ( y - 1 



and the masses 



312 = \ V 7 ^) Qv^+l) 

mi = y/Pi/i and m 2 = -^/W 4 ■ 



(61) 



(62) 



We notice that the mass mi is positive, while the mass m 2 is negative. Introducing the momenta tt p i = —id/ds p i 
and 7T P 2 = —id/ds P 2, we eventually obtain the Hermitean Hamiltonian 



fcl 2 «2 2 

E >1 + 7T p2 
2mi ^— ^ 9<m„ 



911 



with now canonical conjugate variables, [s 9 i,7r p:/ ] = iS pq 6ji. In second quantized form it reads 

J ZTTi^ , . J \X X ) 



322 



(Spl - S«jl) 2 < (S p2 - S g2 ) 2 



£ 



312 



(Spl 



(63) 



(64) 



The Hamiltonian l|63|) describes a one-dimensional interacting many-body system for two kinds of k\ particles at 
positions s p i, p — 1, . . . ,k\ and & 2 particles at positions s p2 , p — 1, . . . , fc 2 on the s axis. 

The supcrunitary model in the form 1631) may be employed to describe electrons in a quasi-one-dimensional semi- 
conductor, see Fig. |21 The electrons are subject to a periodic potential. There is an upper and a lower band, separated 
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FIG. 2: Electrons in the upper (black circles) and lower (open circles) band of a quasi-one-dimensional semiconductor. The 
dispersion relations as function of the wave number k are indicated by the parabola and the inverted parabola. The particles 
are then mapped onto the s axis in the bottom part. 

by a gap. The electrons in the upper band have a positive (effective) mass, while the electrons in the lower band 
close to the gap have a negative (effective) mass. This is due to the dispersion relation e& as function of the wave 
number k. Its second derivative, i.e. the inverse mass, is positive in the upper, but negative in the lower band |49j|. We 
recall that the coupling constants <?y are not arbitrary, they are functions of both (3\ and 02 ■ This makes it possible 
to model repulsive as well as attractive interactions between equal particles and also between different particles by 
choosing proper parameters (3\ and /3 2 . We mention that the spectrum has to be bound from below by an additional 
mechanism if one wants to derive thermodynamical quantities. 



B. Orthosymplectic Model 



As the orthosymplectic model <|34[) is derived from the symmetric superspace OSp (fci/2fc 2 ), it has additional 
symmetries, comprising the ones found in the models based on the ordinary groups O(JV) and Sp(27V). There is 
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a symmetry of point reflections about si = and about S2 — 0. This renders the differential operator of the 
orthosymplectic model (|34|l real and thus Hermitean as it stands. It describes a quasi-two-dimensional physical 
system. One set of particles at positions s p i is confined to the s± axis and a second set of particles at positions s P 2 
confined to the orthogonal S2 axis. As in the superunitary model, all particles interact through a distance dependent, 
inverse quadratic potential. The point reflection symmetry about the two axes implies that each particle at the 
position s p j > with the momentum ir P j has a counterpart at the position — an with the momentum —ir p j- Moreover, 
due to the reflection symmetry, the particles are also subjected to a confining or deconfining inverse quadratic central 
potential. This generalizes the situation described by the models from the ordinary groups O(N) and Sp(JV) p^ |. 

However, the orthosymplectic model has yet another important feature. Closer inspection reveals that the potentials 
also contain angular dependent terms. We now show that these are dipole-dipole interactions, referred to as tensor 
forces in nuclear physics |50j| . The general form of such a dipole-dipole interaction in d dimensions reads 



V(f p ,f q ) = v(\f p - r q \) (e pq ■ <7 p )(e 



-VI 



(65) 



where r p is the position of particle p and S p the dipole vector attached to it. The vector e pq is the unit vector pointing 
in the direction r p — r q . The potential v(r) depends on the distance between the particles only. In nuclear physics, it 
is short-ranged |50j . in our case the potential comes out inverse quadratic, v(r) = 1/r 2 . In the following discussion 
we assume d — 2. This assumption is not a necessary one. Interpretations in higher dimensions are also possible, 
but may be discussed elsewhere. We notice that the functional form of the potential, when derived from a Poisson 
equation, depends on the number of spatial dimensions. Thus, one should not view the dipole-dipole interaction as 
stemming from a Coulomb potential in the present two-dimensional interpretation. For d = 2 we write Eq. Ij65|) more 
explicitly as 



GpCJq 



-e T 

|2 PQ 



cos(i? p + d q ) 
sin z9„ cos i?„ 



cos 'dp sin i9 q 

- COs(#p + 'dq) 



(66) 



with dp = (7 p (cos "dp, sin^p). In our quasi-two-dimensional model, there are three possibilities for the distance vectors 
r p — r q . Expressed in the coordinates s p \ and s P 2, they read 



±S pl ± S ql 









±Spi 







_±S p2 ± S q2 _ 
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_±S q 2_ 



(67) 



depending on which axis the particles p and q move. The angular dependent interaction in Eq. I|34fl can easily be cast 
into the form (|66|l . 

Hence, the orthosymplectic model (|34[l describes the motion of two kinds of charged particles with dipole vectors 
attached to them. The interaction comprises, first, a central potential, second, an only distance dependent potential 
and third a tensor force. Two examples are sketched in Fig-EI Restricting ourselves to even k\, we cast the Hamiltonian 
into the new form 



kl 2 2k 2 2 

2mi ^ 2m2 ^ ( Spl 



P =i 



E 



h-12 



2m 2 



/in 



4i + s 2 



<f2 



p=l *P2 



= 1 S P1 



'22 



P~^ (Sp2 - S q2 ) 2 



E 



(e pq ■ (7i) (e pq ■ (72 ) - (?i • <?2/2 



(68) 



5 pi 



5 g2 



and match it on Eq. Ij34|) by adjusting the parameters. The masses are uniquely determined. They are now both 
positive and given by 



mi = VPi/4 and m 2 = 



(69) 



In order to determine the other free parameters in Eq. I|68(l we have to choose specific directions of the dipoles. 
There are various constraints. All dipoles attached to the particles on the negative si axis must point into the same 
direction, described by the angle i?i_, say. Similar constraints apply to the dipoles on the other half-axes. We denote 
the corresponding angles by $1+ for the positive si axis and with $2- and $2+ for the half-axes in S2 direction. 
Nevertheless, the four angles can not be chosen arbitrarily, there are some further constraints which are given mlCI 
together with a complete list of all possible combinations of different directions. Here we only consider the possibility 
■i9i_ = $ 1+ = $1 and i? 2 - = $2+ = i9 2 - Moreover, there is some arbitrariness for choosing the moduli <7j, j — 1,2. 
For the sake of simplicity, we assume the strengths of both dipoles to be the same <j\ = 02 — 
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J 



FIG. 3: Two realizations of the orthosymplectic model. Left: 2fei = 4 particles on the si axis and 2/c2 = 6 particles on the S2 
axis. The dipole vectors on the same axis have the same direction. Tensor forces are indicated as thicker and thinner dashed 
lines, corresponding to the strength of the force. The central forces and the distance dependent forces are not depicted. Right: 
a case with different directions of the dipole vectors on different sides of the same axis. 



The strength of the central potential in the Hamiltonian (|68|l is given by 



J 8 V 2 



(70) 



When trying to determine the coupling constants h%j in the Hamiltonian JBHJl, we face yet another type of arbitrariness. 
There are at least two possibilities. The tensor force could, first, act between pairs of particles one on either axis 
or it could, second, acts between all particles. We choose the second option as it seems more natural. The coupling 
constants are then given by 



h n = 
h 2 2 = 
h\2 = 



(T 2 COS 2??i 

a 2 cos 2$ 2 



(71) 



(72) 



The strength of the dipoles is determined through the relation 

A sketch of two possible realizations is given in Fig. [3J Notice that the tensor force between two dipoles vanishes at 
a relative angle of 45° between the particle positions. 

Of course H in Eq. l|fjg|) and the operator H, say, on the left hand side of Eq. are still not identical. For H 

and H to be equivalent, the time evolution for the many-body wavefunction ipj^ 1 ^ 2 (s, t) has to be the same. Thus, 
the corresponding time dependent Schrodinger equations have to fulfill 



d 



•_ ,(/3i,/3 2 ) f t) 



(73) 



Thus, the wave function at t — must already have the reflection symmetry 

<;fef 2) ( S ,0)=<^f ) (- S ,0) at t = 0. (74) 

The different interaction strengths are sketched by different widths of the interaction lines. In [21 all possible combi- 
nations of the dipole directions are derived. They are shown in Fig. Apart from some sign changes, all formulae 
given above for the coupling constants are valid for odd ki as well. 
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FIG. 4: Panel of the possible dipole directions in the orthosymplectic model as derived Icl 



VI. SUMMARY AND CONCLUSIONS 

Using supersymmetry, we derived new classes of models for interacting particles. We obtained, first, a superunitary 
model which is based on the supergroup GL(fci/fc2) and on the symmetric superspace GL(fci/2fc2)/OSp(fci/2fc2) and, 
second, an orthosymplectic model which is based on the supergroup OSp(fci/2/c2). It is crucial that these models 
depend in a non-trivial way on two real parameters (3\ and 0%. Our models extend and include the models of the 
CMS type in ordinary space. 

Moreover, our superunitary model contains the supersymmetric constructions derived in Refs. |3Ct l3l|. The latter 
depend on one parameter o nly, implying that they are defined on a one-parameter subspace in the two-dimensional 
{01,02) plane. In Refs. [32l l33l]. an ad hoc construction of models for different kinds of particles was given, no con- 
nection to supersymmetry was established. Not surprisingly, our superunitary model is recovered for some parameter 
values in this construction. In our approach, the connection to supersymmetry is the essential point. It allowed us 
to explicity construct a complete set of solutions in terms of recursion formulae for a trivial and a non-trivial one- 
parameter subspace in the (0i,02) plane. This strongly corroborates the hypothesis of exact integrability. However 
an ultima te p roof is still lacking. The non-trivial one-parameter subspace coincides with the space considered in 
Refs. [3(3, l3l| . In these studies, solutions in terms of deformed Jack polynomials were derived. The relation of the 
recursion formula derived here and the deformed Jack polynomials has to be further investigated. The recursion 
formulae seem to be generating functions or, equivalently, proper resummations of the deformed Jack polynomials. 
Recursion formulae on other one-parameter subspaces are likely to exists. It would be most interesting to gain deeper 
insight into the role of the one-parameter space where solutions have been worked out. Work is in progress. 

We showed that our models have a very natural interpretation. The superunitary model describes electrons in the 
upper and lower band close to the gap in a quasi-one-dimcnsional semiconductor. The orthosymplectic model applies 
to a quasi-two-dimensional system of two kinds of particles confined to two orthogonal directions. Dipole vectors are 
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attached to the particles. The interaction consists of central, distance dependent and tensor forces. 
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APPENDIX A: CALCULATION OF THE COMMUTATION RELATION l|45|) 

Since the parameters Si 2 and the integration variables s' i2 are related through the linear relation of Eq. I|38[) we have 
the differentiation rules for an arbitrary function /. 

<} /(4) = Jrf(4), i/W = o i = i...fe. (Ai) 



ds l2 ds' i2 ' ds i2 



Acting with onto the integral yields 



The last term in the integral has to be integrated by parts using the rule 



Ig%|2 ^-/(4)|(A2) 



(is p2 - is q2 ) 3 ds' q2 



d[m P \ 2 J^rf(4) = I f(4)d[Z P ] ■ (A3) 



>2 



We obtain 

Ag» / Wi&tW = / 45]#A<f /«) - 



W „(/3) 



-(A4) 



The proof is complete if the second integral vanishes identically. It is a straightforward exercise using the definition 
of in Eq. H43|) and identities such as 



ISgl Iggl IsfcJ = 



( _ , _. fc (is P 2 - is q2 ) 3 (is p2 - is k2 ) 2 (is q2 - is p2 ) 2 (is q2 - is k2 ) 

= (A5) 

to show that this is so. 



APPENDIX B: DERIVATION OF THE PROPERTIES (H71) 

We restrict ourselves to the proof of the second equality Eq. I|47(l . The proof of the first one is along the same lines 
but much simpler. It is useful to introduce the operators 



A (c« - —B {c ' 0) (s)— (Bl) 
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such that 

A M) + A (c,/3) = A (c,/3) _ (B2) 

In Eq. (IE} B^(s) is the function B [ ^ h) (s) of Eq. J23 on the hyperbola f3i = 4/f3 2 . It also proves useful to 
split up the the "mixed" part of the measure function H%p\s, s') as follows 

VBF ( S = S = / i BFl( s l J s 2,s' 2 )MBF / i ) (si,S2,S2) 

= / i BF2 ) ( s l J s i. s 2)MBF2 ) ( s l' S l: s 2) • ( B 3) 



For reasons of clarity we suppress in the following the arguments of the measure functions. Hitting the integral with 

;et 

a< c /> / d^j(s') = 



Aj F we get 



A (c,/3) (i (c,/3) 
a sF r l BF2 



(c/j)2i //3 



E 



1 



«Sj2 - Sfci iSj2 



^r)^*^/*^"^ (B4) 



Here and in the sequel we use the convention that operators in squared brackets act only onto the functions inside 
the squared brackets. We now pull the differential operators in s P 2, p = 1 . . . &2 into the second integral using the 
identity Eq. (|45|l and the differentiation rules Eq. 1|A1|I . We obtain 



Ai^ Jd^f(s') = jd^»([fiP (A^ + Aif)^ 

2i f/3 



(c,/3) 



[3 \2 



1 S 

2« f/3 
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/3 V2 



lS j2 ~ S k l lS j2 ~ S kl 
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(9 s 



J2 



(c,/3) (c,/3) 
"Mf ' ^BF 



/ A*bf 

J 2 



- s k i iSj2 - 



(B5) 



j y--j* --j* °kl/ j2 

The last term in the right hand side has to be integrated by parts using the rule Eq. (|A3(I . Then we can write 



A 



sF 



(c,/3) 
M_BF 



^1 5 



(c,/3) 



J"2 



(B6) 



where Mp(s, s') is a rather lengthy expression, which contains no further derivatives. In order to yield the calculations 
traceable we state the expression explicitly 
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For A^ 3 - 1 we proceed analogously. In this case we use the identity Eq. l(4"5|) for /is(sx,sV) and Ai^- 1 respectively, 
which has been proved in Ref. 39]. We also need the integration formula (4.7) of Ref. '27\. The outcome can be 
written in the same form as Eq. (|B6|I 



A (c,/3) 



^BF a / ^BF 



1 2 



/4t )_1 (Afe« + Ai?f ) + M fl ( Sl ✓)) f(s') 



Here Mg(s, s') is again a rather unhandy expression 
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^ V(**ia-«fci) a (*Si2-4i) 2 
Adding Eqs. (|B6|) and (|B8|> yields the desired result, provided that 

M B (s, s>) + M F (s, s>) = [/4t rl + 



fci-i 



„(c,/3)- 1 ^ „(c,/3) 
OS jl 



—y 



(c/3)" 1 9 ( C)/3 ) 
^BF a j PBP 
CSj2 



(B8) 



(B9) 



(BIO) 



With the definitions of Mb, Mp in Eqs. QB7fl. (|B9|I and of in Eq. I|43|) it is a straightforward but extremely 

tedious exercise to show that Eq. i|B10|l is true. 



APPENDIX C: CONSTRAINTS AND POSSIBLE CHOICES FOR THE DIRECTIONS OF THE DIPOLES 

IN THE ORTHOSYMPLECTIC MODEL 

In order to match the right hand side of Eq. I|34|) with the Hamiltonian Eq. I|68|) the four angles $i±, $2± of the 
dipoles have to meet the following three conditions 

= cos(2i?i_) + cos(2tf l+ ) - 2 cos(iV + , i = 1, 2 

= sin(i? 1+ + a _) + sin(^i_ + tf 2+ ) - sin(i?i_ + tf 2 _) - sin(tf 1+ + § 2+ ) . (CI) 

The first two equations have three solutions each. One solution is di+ — The other solutions are at = \ — 
and at = ^ — i = 1,2. Whenever $i+ = is chosen as solution the third equation in Eq. (|C1|I does not 
yield any further condition. Then two of the four angles can be chosen arbitrarily. In general either one or two angles 
can be chosen freely depending on which solution is selected. All possibilities are compiled in Table 1. 

The two cases depicted in Fig.Olcorrespond to the entries (1, 1) and (2, 3) in Table 1. In Fig.^Ja typical configuration 
for each entry of Table 1 is depicted. We also state the general formulae for the coupling constants hij and fi restricting 
ourselves to fci even, 

hu = (y- 1 ) + ^(cos(2tfi + ) + cos(2i?i_) + 2cos(2tf 1 _ + tf 1+ )) 

fi = (-l)' 1 ^ (y- 1 ) + j| ( cos ( 2 ^i+) + cos(2^_) + 2 cos(2^_ + 
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TABLE I: The angles which can be chosen freely depending on the selected solution of Eq. 1CH , 





■&2+ = 02- 


0a+ = f - 02- 


02+ = ^ - 02- 


01+ = 01- 


01-, 02- 


01-, 02- 


01-, 02- 


01+ = f - 01- 


01-, 02- 


0i_ , (0 2 - = f,f ) 

#2- , (01- = f,f) 


01-, 02- 


01+ = ^-0i_ 


01-, 02- 


01-, 02- 


02- , (^_ = ^,2f) 



In the general case, we find 

(Ti(T2 (COS(0 1+ + 2+ ) + COS(0 1+ + 2 -) + CO S (0!_ + 2 +) + CO S (0!_ + 2 -)) = (VPl ~ fa) (VPiA + (C3) 

for the moduli squared of the dipoles. 
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